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(54) Title: METHOD AND APPARATUS FOR MEASURING SHAPE OF OBJECTS 
(57) Abstract 

There is disclosed a method for mea- 
suring the shape of objects which might have 
surface discontinuities using projected fringes, 
m which method the pitch of the fringes is 
varied over time and a sequence of phase- 
stepped images is recorded from which a three- 
dimensional complex intensity is calculated 
and either the phase of the distribution is un- 
wrapped along the time axis or the Fourier 
transform of the distribution is calculated alone 
the time axis. 6 
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METHOD AND APPARATUS FOR MEASURING SHAPE OF OBJECTS 
This invention relates to measuring the shape of objects. 

The use of projected fringes for measurement of surface profile is a well- 
developed technique as described for example by V. Srinivasan, H.C. Liu and M. 
Halioua, "Automated phase-measuring profilometry of 3-D diffuse objects » Appl Opt 
23 3105-3108 (1984). Parallel fringes are projected onto the object surface, either by a 
conventional imaging system or by coherent light interference patterns. If the projection 
and observation directions are different, the phase distribution of the measured fringe 
pattern includes information on the surface height profile of the object. The sensitivity 
can be improved by viewing through a line grating: the technique is then referred to as 
projection moire. 

Automated analysis of the fringe patterns is normally carried out either by 
the Fourier transform method (M Takeda, H Ina and S Kobayashi, "Fourier-transform 
method of fringe-pattern analysis for computer-based topography and interferometry" 
J Opt. Soc. Am 72 156-160 (1982)), or by phase-stepping of the fringe patterns (J H 
Brunmg, D.R. Herriott, J.E. Gallagher, D.P. Rosenfeld, A.D. White and D.J. Brangaccio 
"Dignal wavefront measuring interferometer for testing optical surfaces and lenses" 
Appl. Opt. 13 2693-2703 (1974)). Both produce wrapped phase maps (i.e., phase values 
lymg ,n the range -« to +n); the 2n phase jumps must then be removed by the process 
known as phase unwrapping to recover the surface shape. Phase unwrapping is normallv 
carried out by comparing the phase at neighbouring pixels and adding or subtracting 
multiples of 2. to bring the relative phase between the two pixels into the ranee -„ to +, 
This causes problems when the technique is applied to real engineerine objects since 
such objects often contain sudden jumps in profile. It can then become impossible to 
unwrap correctly across such real discontinuities and large (multiples of 2.) phase errors 
can propagate across the image. The problem can be reduced in principle by measuring 
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several phase maps at different sensitivities. Only one unwrapped map will then be 
consistent with all the wrapped maps. White light interferometry can also measure 
absolute surface heights, but is difficult to apply when the height variations exceed the 
range of piezo-electric translators (PZT), typically a fraction of a mm. 

The present invention provides methods and apparatus which can rapidly 
and unambiguously measure the profiles of discontinuous objects. 

The invention comprises a method for measuring the shape of objects which 
might have surface discontinuities using projected fringes, in which the pitch of the 
fringes is varied over time. Each time the fringe pitch is changed a sequence of phase- 
stepped images of the fringe pattern is recorded, from which a complex intensity map is 
calculated. The sequence of such maps forms a three-dimensional complex intensity 
distribution. The phase of this distribution is then unwrapped along the time axis or else 
the Fourier transform along the time axis is calculated to yield data representative of the 
shape of the surface. 

Intensity maps may be recorded as pixel images. 

The pitch of the fringes may be varied by piezo-electric device means. 

The fringes may be produced by Michelson interferometer means, and the 
phase and spatial frequency of the fringes may then be varied by plane mirrors forming 
part of the interferometer means. 

The fringes may, however, be generated by interference between laser light 
traversing dissimilar paths of which at least one is variable for changing phase and spatial 
frequency. 
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The fringes may also be generated from incoherent light by a spatial light 
modulator (SLM) which can be operated to change the phase and spatial frequency of the 
fringes. 



The invention also comprises apparatus for measuring shape of objects 
which might have surface discontinuities using projected fringes comprising fringe 
generating means projecting fringes on the object to be measured, phase stepping means 
to vary the phase of the fringes and fringe pitch varying means adapted to vary the pitch 
of the generated fringes over time, recording means adapted to record the sequence of 
phase-stepped images ,as the fringe pitch is varied over time and computing means 
programmed to calculate a three-dimensional complex intensity distribution either to 
unwrap the phase of this distribution along the time axis or to calculate its Fourier 
transform along the time axis to yield data representative of the shape of the surface. 

The change in the fringe pitch between any two successive phase maps must 
be sufficiently small that the resulting unwrapped phase change lies in the range -n to +. 
for all points within the measurement region. This is equivalent to the requirement that, 
provided the position of the central fringe remains unchanged during the changing of the' 
fringe pitch, no more than one fringe is added to the projected fringe pattern between the 
acquisition of two successive phase maps. The requirement may however be relaxed in 
certain circumstances. 

The fringe generating means may comprise Michelson interferometer means 
in which the fringe phase and pitch varying means comprise plane mirror means. The 
interferometer means may comprise a beamsplitter cube with orthogonal adjacent faces 
and adjustable mirrors in juxtaposition to said adjacent faces, of which one mirror mav 
be adjustable towards and away from its associated face (for phase chaneine) and the 
other adjustable as to tilt (for spatial frequency changing). Alternatively the tilt and 
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translation operations may be combined into a single mirror mount. At least one mirror 
may be adjustable by means of a piezo-electric device. 

The fringe generating means may comprise laser illuminated optical fibre 
means providing two dissimilar light paths. The phase varying means may comprise 
optical fibre length varying means, which may comprise a piezo-electric cylinder around 
which the fibre is wrapped and which is electrically actuated to vary its diameter. A 
second piezo-electric device can relatively shift adjacent light output ends of the optical 
fibres to change the fringe spatial frequency. 

The fringe generating means may comprise an incoherent, white light 
source and a spatial light modulator driven by an external signal source so as to produce 
synthetic fringes such that the transmitted light intensity has an approximately sinusoidal 
profile. The arrangement may comprise a condenser lens system and a projection lens 
projecting the fringes onto the object to be measured. 

The transmittance of the spatial light modulator may be wavelength- 
dependent thereby allowing several fringe patterns of differing phase or spatial frequency 
to be projected simultaneously as patterns of differing colour within a single projected 
image. 

The recording means may comprise a pixel image recorder such as a video 
system in which the object is imaged by a lens on to the light-sensitive image detector. 
The sensitivity of the image detector may be wavelength-dependent thereby allowing 
several fringe patterns of differing phase or spatial frequency to be recorded as a single 
coloured image. 

The computing means may comprise a processor arranged to operate in real 
time by carrying out simple arithmetic operations on an entire pixel image frame 
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simultaneously. Parallel processors or a fast serial processor could be used. It is of 
course not always required to operate in real time. The processor may have sufficient 
random access memory to store the entire three-dimensional complex intensity 
distribution so that the one-dimensional Fourier transform can be calculated along the 
time axis. 



have discontinuities using projected fringes in accordance with the invention will now 
be described with reference to the accompanying drawings, in which 



Methods and apparatus for measuring the shape of objects which might 



Figure 1 



is a wrapped phase map of an object (grey shades ranging 
from black to white represent phase values ranging from -u 



to +7i ; 



Figure 2 



is the result of unwrapping the phase map of Figure 1 
through a sequence of incremental phase maps; 



Figure 3 



is a diagrammatic illustration of a first measurement 



apparatus; 



Figure 4 



is a diagrammatic illustration of a fringe generating means; 



Figure 5 



is a diagrammatic illustration of another fringe generating 



means; 



Figure 6 



is a diagrammatic illustration of a second measurement 
apparatus; 
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is a representation of the measured intensity at three pixels 

(a) , (b) and (c) as a function of phase shift index (prvertical 
axis) and time (t:horizontal axis) in the measurement volume 
of Figure 6; 

Figure 8 shows the unwrapped time histories for the three pixels (a), 

(b) , (c), of Figure 7; 

Figure 9 comprises contour maps of the 2-D Fourier transforms of the 
intensity distributions for the three pixels shown in Figure 7, 
the horizontal and vertical axes being k, and k p9 representing 
spatial frequencies along the time (t) and phase shift index 
(p) axes; and 

Figure 10 shows cross-sections of the contour maps of Figure 9 along 
the line k p = 1 - the vertical axis being 500 arbitrary units 
corresponding to the contour interval in Figure 9. 

Figure 1 is a typical wrapped phase map such as would be obtained by 
conventional projected fringe measurement of a test object consisting of three 4mm wide 
vertical grooves on 6mm centres machined in a flat piece of plastic. The edges of the 
first and third grooves are visible on the left and right sides of the image. The depths of 
the three grooves are approximately 1mm, 2mm and 4mm. It is impossible to unwrap 
this phase map correctly by conventional spatial methods since the phase jumps at the 
groove edges fall outside the range (-it, +*). The relative depths of the grooves are 
therefore uncertain to integral multiples of 2n. The test object was illuminated along the 
grooves to avoid shadowing. 
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The problem is resolved by applying the technique of temporal phase 
unwrapping ( J M Huntley and H O Saldner, Temporal phase-unwrapping algorithm for 
automated interferogram analysis, Appl. Opt. 22, 3047-3052 (1993)). However, ways 
must be provided to vary both phase and spatial frequency or pitch of the fringes. 

Figure 3 illustrates one such way, using a Michelson interferometer 
arrangement. Coherent light from the laser is directed via spatial filter SF and objective 
lens OL, to a beamsplitter cube BS whence fringes are directed at the measurement 
region R. The fringes are imaged in a CCD camera via objective lens OL 2 . 

M, and M 2 are plane mirrors moved by piezo-electric translation devices 
PZT, and PZT 2 . M 2 is juxtaposed to one face of the cube and moved towards and away 
therefrom by its PZT to alter the phase of the fringes at R. M, is juxtaposed to an 
adjacent, orthogonal face of the cube and tilted by its PZT about fulcrum F to alter the 
spatial frequency of the fringes. 

In the method, M, is tilted to adjust the angle between the two beams in a 
number of steps, for example, 22 steps. For each angular step, four phase-stepped 
images are required. A phase change of n/2 is introduced between each of the successive 
four images by moving mirror M 2 . 

Using an 80mW frequency - doubled diode-pumped YAG laser (A. = 
532nm) with a 100mm focal length for objective OL, at f73.5, a speckle diameter of 
3.5um is obtained. The active pixel area of the CCD camera is 5.8 x 7.4um- giving 
several speckles per pixel, which substantially reduces the effect of speckle noise. 

Angle a (Figure 3) was set to 20» for making the image of Figure 1. 
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The succession of images obtained is transferred to a fast computer for 
analysis. Incremental phase maps are determined using Equation 9 from Huntley and 
Saldner, loc cit (Equation (3) below). These are added by Equation 10 (Equation (1) 
below) to give the unwrapped phase map. After subtraction of a uniform phase ramp, 
the unwrapped phase map is as shown in Figure 2. The surface profile has been 
measured accurately and unambiguously. 

Figure 4 shows an arrangement in which laser light is directed along a fibre 
41 to a coupler so that coherent, in-phase light travels with equal amplitudes in fibres 
41 A, 4 IB to recombine at their ends to produce interference fringes. Fibre 41 A is 
wrapped around a piezo-electric cylinder PZT, which expands and contracts radially 
according to an applied voltage, stretching fibre 41 A to different extents, thereby 
shifting the relative phase of the recombining beams. The ends of the fibres 41 A, 4 IB 
are relatively shifted in a direction perpendicular to the fibre axis by a second piezo- 
electric device PZT 2 which changes the fringe spacing. 

Figure 5 shows an arrangement using white, incoherent light from a source 
5 1 passing through a condenser 52 to a spatial light modulator 53 controlled by a signal 
generator 54 and a projection lens 55 to produce variable fringes. 

The various embodiments of apparatus can be built into a comprehensive 
system with self-contained computing to produce various instruments adapted to make 
different kinds of measurement. The instruments will facilitate human body 
measurement for clothes manufacture and even bespoke tailoring, quality and process 
control on production lines, robot vision, medical and dental measurements, e.g. for 
prostheses, duplication of e.g. delicate artefacts, verification that a CAD/CAM product 
matches the design specification, reverse engineering and non-contact interfaces for 
virtual reality systems, among many other applications, 
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The calculation of the phase maps can be done in many different ways. 
Some use only three phase-stepped images, as compared to the four used in Huntley and 
Saldner, loc.cit, but using five or more can lead to greater accuracy. 

Various modifications of the above-described techniques are possible with 

advantages. 

Instead, for example, of projecting a series of fringes in the sequence t - 1, 
2, 3, 4, 5, ...5, (t being the number of fringes across the field of view) a sequence t=l, K, 
K 2 , K 3 , ...s is projected where K is any number (not necessarily integral) greater than 1, 
for example, 2. 

This results in an improvement in the success rate of the unwrapping 
procedure as well as a reduction in the time required for image acquisition and 
computation. With K = 2 there is a good compromise between unwrapping reliability 
and computational cost. 

Intermediate phase values can be incorporated into the calculations by linear 
least squares analysis to reduce the root-mean-square height error. 

By way of illustration, and with reference to Figure 6, the fringe pitch can 
be changed over time so that initially a single fringe is projected from the spatial light 
modulator SLM into a field of view region in the measurement volume so that the fringe 
phase ranges from -« on one side to +n on the other side. For subsequent time values 
(/ = 2, 3, ...s) the number of fringes is also set equal to t so that the phase range increases 
to (- tn , tiz). Figure 6 shows the situation for the case t = 2. For each t value, a number 
of phase-stepped images is acquired [p = 1 , 2. ...q). The measured intensity, / is therefore 
a function of four variables : pixel coordinates {m, n) in the field of view, p and t. 
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Figure 7 shows the / (p, t) map for three pixels with q = 4 and phase steps 
of 7r/2. The first pixel (Figure 7 (a)) receives light from a point in the field of view where 
the phase of the projected fringes is increasing with time. Figure 7 (b) corresponds to a 
plane where the phase is not changing appreciably over time - in the middle of the 
projection field of view, for example, with the fringes coming from either side of it. 
Figure 7 (c) corresponds to a location on the opposite side where the phase is decreasing 
with time t. 

The phase calculations for any point (m 9 n) involve computing the time- 
varying unwrapped phase by summing the wrapped phase changes between successive 
phase maps : 

*(/)=2AO//-l) (1) 

where A*(/ 1 /)=*(0-*(/) (2) 

and fc(0) is defined to be zero. Subscript w denotes a wrapped phased value (i.e. one 
lying in the range (-*, *)), and subscript u denotes an unwrapped phase value (i.e. 
adjusted by addition of the correct integral multiple of 2n). The phase changes Afc w (ij) 
are calculated directly from the two sets of intensity values at times t=i and t-j using the 
formula: 



Atfjy) = tan 1 



A/J0A/ l3 (/>-A/ l3 (i)A/ 42 (/) 



A/ n (0A7 I3 (/-)-A/ 42 (0A/ 42 a) 



(3) 



where 

M u (t)=I(kj)-I(U). (4) 



The total unwrapped phase, Y, given by 
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(5) 



then defines the plane which contains the scattering point. Figure 8 shows the 
unwrapped phase calculated in this way from the intensity data in Figures 7(a) - (c). 

The unique coordinates in space of the scattering point within the 
measurement volume are then defined by the point of intersection between this plane and 
the line passing from the pixel (m. n) through the centre of the camera lens into the 
measurement volume. In practice, the surface of constant phase change may deviate 
from a plane due to imperfections in the fringe generating means. Similarly the line 
passing from the pixel {m, ») through the centre of the camera lens into the measurement 
volume may deviate from straightness due to imperfections in the recording means. Both 
of these errors may be corrected for by calibration of the apparatus. 



most 



The unwrapping can be carried out with fewer phase maps, and under 
conditions more reliably, by taking only the phase values in the sequence / = 1,2,4.. .s, 
corresponding to the case K = 2. The general term in this sequence is / = V, where v = 
0, 1 The wrapped phase changes between successive frames are then unwrapped 

using the recursive equations 



(6) 



A*„(2/,0) = A*^/,,) ♦ A* (r,0). 



(7) 



for / - l,2,4,...j/2, and starting from the phase value A* U (1,0) = *„(i). £/{*„*,} 
represents the unwrapping operator which subtracts an integral multiple of 2n from «, 
such that * , - lies in the range -n to +n , i.e. 



(8) 
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where NINT denotes rounding to the nearest integer. 

The measured phase values will in practice have errors associated with 
them. If the errors are additive and Gaussian, the measured linear sequence of 
unwrapped phase values can be written 

= o>r+e^ (9) 

where is a random variable with a zero-mean Gaussian probability distribution. o> 
is a non-dimensional frequency, representing the rate of change of phase with non- 
dimensional time /, i.e. the slope of the graphs shown in Figure 8. 

In the presence of noise, a better approach to the estimation of u> than 
simply using the final unwrapped phase value is a fitting in the least squares sense of the 
line defined by the equation G = coMo all the experimental data for the pixel analysed. 

Least squares fitting does not significantly improve the estimate of u> when 
using an exponentially growing sequence of fringe frequencies. The reason for this is 
that the phase values used are all clustered around the low-/ end of the fc-f graph. Low 
/ values do not provide reliable estimates of the gradient, and so do not contribute 
significantly to the least-squares estimate for <•>. However, the algorithm can be modified 
by starting at the maximum fringe density (r = s), and reducing the number of fringes by 
1,2.4, 8.. .s/2. This ensures that the measured phase values are now clustered at the high-f 
end of the *-r graph. Equations (6) and (7) are modified to 

A* u (5^V-2f , )=t/{A« w (^-r^-2r'),A^ B (w-r')} (10) 



A* l< (5,J-2f')=A^ lJ (5-r^^-2f•)-A^ l^ (^-r , ). 



(11) 
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for t' = 1,2,4,-. -s/2, where the starting phase value A* u (5,j-I)=a<X> w (j,5-1). 

Both the unwrapping reliability and the computation time for this method 
are identical to those for the forward exponential sequence, yet the measurement 
accuracy is improved by a factor of approximately ^/to£>. It is therefore recommended 
as the most appropriate algorithm when data acquisition and processing times must be 
minimised. 

An alternative method of analysing the data which avoids the need for 
unwrapping is based on Fourier analysis of the 2-dimensional intensity distribution 
measured at each pixel. The 2-D transforms of the intensity distributions from the three 
pixels shown in Figure 7 are presented in Figure 9. The horizontal and vertical axes 
represent spatial frequencies ^ and along the / (time) and p (phase-shift index) 
directions, respectively. The transform consists of three peaks: a dc term centred on k x 
= k p = 0, and two sidelobes symmetrically oriented about the dc peak. If the phase- 
shifting algorithm used covers one complete cycle, then the sidelobes will always appear 
on the lines k p = ±1 . The position of the peak along the line k p = 1 is then proportional 
to the rate of change of phase as the fringes sweep past the pixel of interest, i.e. is 
proportional to cj. Figure 10 shows the cross section k p = 1 for the three transforms of 
Figure 9. 

In practical terms, the curves shown in Figure 10 can be calculated directly 
by forming a 1 -dimensional complex array 



M0=A/ 13 (0+iA/ 42 (r) (12) 
where / is the square root of -1, and then taking a one-dimensional Fourier transform 
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= S Kr ♦ I )exp( -Inikfts} ( 1 3) 

Computationally, the aim is to find the value k x = K which maximises 
(//ft)) 2 . This may be done as a two-step process: firstly, Equation (13) is evaluated for 
integral values of k, using the Fast Fourier Transform algorithm; the value k t which 
maximises \H(k^\ 2 may then be used as the starting point for an optimisation based on 
a bounded Newton-Raphson algorithm, in which k, is no longer constrained to be 
integral. The random error from the method can be reduced further by noting that the 
signal peak, which is in general complex, can be made real by multiplying h(t) by exp(- 
2n/ic/s). A final iteration to maximise {Re[//#,)]} 2 jhen excludes noise contributed by 
the imaginary part of the transform. 



Equation (12) can be generalised for use with phase-stepping algorithms 
other than the four-frame one (q = 4) used here. All such algorithms result in two 
numbers : one proportional to cos(*), which is used to form the real part ofh(t), and the 
other one to sin(fc), which forms the imaginary part oih(t). 

The computation time for this method is o($log 2 .y) which is a factor log,* 
longer than that for the linear fitting method. However, the Fourier transform method 
is interesting for two significant reasons; firstly the transform concentrates the signal 
with respect to the noise. The greater the value of s> the higher the signal to noise ratio 
in the vicinity of the peak. Thus the reliability is high and furthermore increases with 
increasing s, unlike the unwrapping methods for which the reliability decreases as s 
increases. Secondly, multiple reflections can occur if the object being profiled has a 
smooth surface. The light received at a given pixel then comes from two or more 
different points on the surface, having different values of <o. For the phase unwrapping 
methods the result is a ripple on the <&-/ graph which causes systematic errors in the 
gradient estimation, whereas the Fourier transform method will result in two or more 
peaks which do not interfere provided the frequencies are separated by at least k, = ±1 . 
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Furthermore, the direct and indirect reflections could in principle be distinguished by 
making measurements of the object from different viewing or illumination directions and 
checking which of the peaks give consistent position measurements. 
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CLAIMS 

1 * A method for measuring the shape of objects which might have surface 

discontinuities using projected fringes, in which method the pitch of the fringes is varied 
over time and a sequence of phase-stepped images is recorded from which a three- 
dimensional complex intensity distribution is calculated. 

2. A method according to claim 1, in which the phase of the distribution is 
unwrapped along the time axis. 

3. A method according to claim 1, in which the Fourier transform of the 
distribution is calculated along the time axis. 

4. A method according to any one of claims 1 to 3, in which the phase maps 
are recorded as pixel images. 

5. A method according to any one of claims I to 4, in which the pitch of the 
fringes is varied by piezo-electric device means. 

6. A method according to any one of claims 1 to 5, in which the fringes are 
produced by Michelson interferometer means. 

7. A method according to claim 6, in which phase and spatial frequency of the 
fringes are varied by plane mirrors. 

8: A method according to any one of claims 1 to 7, in which the fringes are 

generated by interference between laser light traversing dissimilar paths of which at least 
one is variable for changing phase and spatial frequency. 
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9. A method according to anyone of claims I to 4, in which the fringes are 
generated from incoherent light by a spatial light modulator. 

10. Apparatus for measuring the shape of objects which might have surface 
discontinuities using projected fringes comprising fringe generating means projecting 
fringes on the object to be measured, phase stepping means to vary the phase of the 
fringes and fringe pitch varying means adapted to vary the pitch of the generated fringes 
over time, recording means adapted to record a sequence of phase-stepped images and 
computing means programmed to calculate a three-dimensional complex intensity 
distribution. 

11. Apparatus according to claim 10, in which the computing means are 
programmed to unwrap the phase of the distribution along the time axis. 

12. Apparatus according to claim 10, in which the programming means are 
programmed to calculate the Fourier transform along the time axis of the distribution. 

13- Apparatus according to claim 12, in which the fringe generating means 

comprise Michelson interferometer means. 

14. Apparatus according to claim 13, in which the fringe phase and pitch 
varying means comprise plane mirror means. 

15. Apparatus according to claim 14, in which the interferometer means 
comprise a beamsplitter cube with orthogonal adjacent faces and adjustable mirrors in 
juxtaposition to said adjacent faces. 

16. Apparatus according to claim 1 5. in which one mirror is adjustable towards 
and away from its associated face and at least one of the mirrors is adjustable as to tilt. 
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17. Apparatus according to claim 15 or claim 16, in which at least one mirror 
is adjustable by means of a piezo-electric device. 

18. Apparatus according to claim 10, in which the fringe generating means 
comprise laser illuminated optical fibre means providing two dissimilar light paths. 

19. Apparatus according to claim 18, in which the phase varying means 
comprise optical fibre length varying means. 

20. Apparatus according to claim 19, in which said length varying means 
comprise a piezo-electric cylinder around which the fibre is wrapped and which is 
electrically actuated to vary its diameter. 

2 1 . Apparatus according to claim 10, in which the fringes are generated from 
incoherent light by a spatial light modulator. 

22. Apparatus according to claim 21, in which the fringe generating means 
comprise a white light source, a condenser lens system, a spatial light modulator and a 
projection lens projecting the fringes on to the object to be measured. 

23. Apparatus according to any one of claims 10 to 22, in which the recording 
means comprise a pixel image recorder such as a video system. 

24. Apparatus according to any one of claims 10 to 23, in which the computing 
means comprise a processor arranged to operate in real time by carrying out parallel 
simple arithmetic operations on an entire pixel image frame simultaneously. 




SUBSTITUTE SHEET (RULE 26) 



WO 97/36144 



PCT/GB97/00771 



2/7 




LASER 



SF OU BS 



] 0 0 



M- 



PZT. 



PZT- 



Z3 



COUPLER 



LASER 




PZT 2 




SUBSTITUTE SHEET (RULE 26) 



WO 97/36144 „,~^„„,„v„, 

PCT/GB97/00771 

3/7 




SUBSTITUTE SHEET (RULE 26) 



W ° 97/36144 PCT7GB97/0077, 

4/7 





SUBSTITUTE SHEET (RULE 26) 



W ° 97/36144 PCT/GB97/0077, 

5/7 




ao 
4- 



0)<J> 



SUBSTITUTE SHEET (RULE 26) 



W ° 97/36,44 PCT/GB97/0077, 

6/7 




SUBSTITUTE SHEET (RULE 26) 



WO 97/36144 PCT/GB97/00771 

7/7 




SUBSTITUTE SHEET (RULE 26) 



INTERNATIONAL SEARCH REPORT 



Inter. rial Application No 

PCT/6B 97/GG771 



A. CLASSIFICATION OF SUBJECT MATTER. 

IPC 6 G01B11/24 



According to International Patent Classification (IPC) or to both national classification and IPC 



B. FIELDS SEARCHED 


Minimum documentation searched (classification system followed by da. 

IPC 6 G01B 


BsficaQon symbols) 




Documentation searched other than minimum document* Don to the exten 


t that such document* are in< 


:tudcd in the fields searched 



Electronic data base consulted during the international search (name of data base and, where practical, search terms used) 



C. DOCUMENTS CONSIDERED TO BE RELEVANT 



Category " Citation of document, with indication, where appropriate, of the relevant passages 



Relevant to claim No. 



DATABASE INSPEC 

INSTITUTE OF ELECTRICAL ENGINEERS, 
STEVENAGE, GB 

Inspec No. 5556642, February 1997 
SALDNER H 0 ET AL: "Profi lometry using 
temporal phase unwrapping and a spatial 
light modulator-based fringe projector* 
XP082G33884 
see abstract 

& OPTICAL ENGINEERING, FEB. 1997, SPIE, 
USA, 

vol. 36, no. 2, ISSN 0091-3286, 
pages 610-615, 



1,10 



m 



Further documents are listed in (he continuation of box C. 



Patent family members are listed in annex. 



* Special categories of cited documents : 

* A* document defining the general state of the art which is not 

considered to be of particular relevance 
*E" earlier document but published on or after the international 
filing date 

"L* document which may throw doubts on priority daimfs) or 
which is a ud to establish the publication date of another 
citation or other special reason (as specified) 

"O" document referring to an oral disclosure, use, exhibition or 



later document published after the international filing date 
or priority date and not in conflict with the application but 
cited to understand the principle or theory underlying the 



"P* .document published prior to the international filing date but 
later than the pnonty date claimed 



*X* document of particular relevance; the d aimed invention 
cannot be considered novd or cannot be considered to 
involve an inventive step when the document is taken alone 

*Y' document of particular relevance; the d aimed invention 
cannot be considered to involve an inventive step when the 
document is combined with one or more other such docu- 
ments, such combination being obvious to a person skilled 
in the art. 

document member of the same patent family 



Date of the actual completion of the international search 

26 June 1997 



Date of mailing of the international search report 

14.07.97 



Name and mailing address of the ISA 

European Patent Office, P.B. 58 1 8 Patentlaan 2 
NL - 2280 HV Rijswijk 
Td. 31-70) 140-2040, Tx. 31 651 epo nj, 
Fax 31-70) 340- 301 6 



Authorized officer 



Ganci, P 



i PCT/ISA'310 



*h«t) (July 1992) 



page 1 of 2 



INTERNATIONAL SEARCH REPORT 



C.(ConBnujiOon) DOCUMENTS CONSIDERED TO BE RELEVANT 



tnttrr. al Application No 

PCT/GB 97/00771 



Category " Citation of document, with indication, where Appropriate, of the relevant passages 



Relevant to claim No. 



P.A 



DATABASE INSPEC 

INSTITUTE OF ELECTRICAL ENGINEERS, 
STEVENAGE, GB 
Inspec No. 5088483, 

BURTON 0 R ET AL: "The use of carrier 
frequency shifting for the elimination of 
phase discontinuities in Fourier transform 
prof i lometry" 
XP092033870 
see abstract 

& OPTICS AND LASERS IN ENGINEERING, 1995, 
UK, 

vol. 23, no. 4, ISSN 0143-8166, 
pages 245-257, 

APPLIED OPTICS, 

vol. 33, no. 14, 10 May 1994, 

pages 2939-2948, XP000442254 

BURTON D R ET AL: "MULTICHANNEL FOURIER 

FRINGE ANALYSIS AS AN AID TO AUTOMATIC 

PHASE UNWRAPPING" 

see the whole document 

APPLIED OPTICS, 

vol. 34, no. 14, 10 May* 1995, 

pages 2560-2563, XP000511454 

QUIROGA J A ET AL: "PHASE-UNWRAPPING 

ALGORITHM BASED ON AN ADAPTIVE CRITERION" 

see the whole document 



1,10 



1.10 



1,10 



US 5 321 497 A (AI CHIAYU 
1994 

see abstract 



ET AL) 14 June 



DATABASE INSPEC 

INSTITUTE OF ELECTRICAL ENGINEERS. 
STEVENAGE, GB 
Inspec No. 5592018. 

SALDNER H 0 ET AL: "Shape measurement of 
discontinuous objects using projected 
fringes and temporal phase unwrapping" 
XP002033877 
see abstract 

& PROCEEDINGS. INTERNATIONAL CONFERENCE ON 
RECENT ADVANCES IN 3-D DIGITAL IMAGING AND 
MODELING (CAT. NO. 9 7TB 100 134 ) , 
PROCEEDINGS. INTERNATIONAL CONFERENCE ON 
RECENT ADVANCES IN 3-D DIGITAL IMAGING AND 
MODELING (CAT. N0.97TB100134) , OTTAWA, 
ONT. , CANADA, ISBN 0-8186-7943-3, 1997, 
LOS ALAMITOS, CA, USA, IEEE COMPUT. SOC. 
PRESS, USA, 
pages 44-50, 



1,10 



1.10 



Form PCT'lS A/210 (continuation of t*aod ihit() (July 1992) 



page 2 of 2 



INTERNATIONAL SEARCH REPORT 

Lnfomvitjoci on patent family members 



Intei jnal Application No 

PCT/GB 97/00771 



Patent document 
cited in search report 



Publication 
date' 



Patent family 
member(s) 



US 5321497 A 



Publication 
date 



14-06-94 



NONE 





Form PCT/ISA/210 (p*um (unity utnmx) (July 1993} 



This Page is Inserted by IFW Indexing and Scanning 
Operations and is not part of the Official Record 

BEST AVAILABLE IMAGES 

Defective images within this document are accurate representations of the original 
documents submitted by the applicant. 

Defects in the images include but are not limited to the items checked: 

BLACK BORDERS 

IMAGE CUT OFF AT TOP, BOTTOM OR SIDES 
$ FADED TEXT OR DRAWING 

□ BLURRED OR ILLEGIBLE TEXT OR DRAWING 

□ SKEWED/SLANTED IMAGES 

l^jf COLOR OR BLACK AND WHITE PHOTOGRAPHS 

□ GRAY SCALE DOCUMENTS 

□ LINES OR MARKS ON ORIGINAL DOCUMENT 

□ REFERENCE(S) OR EXHIBIT(S) SUBMITTED ARE POOR QUALITY 

□ OTHER: ; . 

IMAGES ARE BEST AVAILABLE COPY. 
As rescanning these documents will not correct the image 
problems checked, please do not report these problems to 
the IFW Image Problem Mailbox. 



